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Quantum-enhanced metrology can be achieved by entangling a probe with an auxiliary system, 
passing the probe through an interferometer, and subsequently making measurements on both the 
probe and auxiliary system. Conceptually, this corresponds to performing metrology with the pu¬ 
rification of a (mixed) probe state. We demonstrate via the quantum Fisher information how to 
design mixed states whose purifications are an excellent metrological resource. In particular, we 
give examples of mixed states with purifications that allow (near) Heisenberg-limited metrology and 
provide examples of entangling Hamiltonians that can generate these states. Finally, we present 
the optimal measurement and parameter-estimation procedure required to realize these sensitivities 
(i.e., that saturate the quantum Cramer-Rao bound). Since pure states of comparable metrological 
usefulness are typically challenging to generate, it may prove easier to use this approach of entan¬ 
glement and measurement of an auxiliary system. An example where this may be the case is atom 
interferometry, where entanglement with optical systems is potentially easier to engineer than the 
atomic interactions required to produce nonclassical atomic states. 


I. INTRODUCTION 


There is currently great interest in quantum metrology: 
the science of estimating a classical parameter <f> with a 
quantum probe at a higher precision than is possible with 
a classical probe of identical particle flux. Given a fixed 
number of particles, TV, the ultimate limit to the sensitiv¬ 
ity is the Heisenberg limit A<f> = 1/N [lj[2]. Naively, the 
choice of probe state is a solved problem; for instance, 
symmetric Dicke states mm and spin-cat states 00 
input into a Mach-Zehnder (MZ) interferometer yield 
sensitivities of \/2/N and 1 /TV, respectively. However, 
in practice achieving quantum-enhanced sensitivities is a 
significant challenge. This is due to both the deleteri¬ 
ous effect of losses [6] and the challenges associated with 
preparing nonclassical states with an appreciable num¬ 
ber of particles jMH- For example, protocols for gen¬ 
erating a spin-cat state commonly require a large Kerr 
nonlinearity, which is either unavailable (e.g. in opti¬ 
cal systems [12]), difficult to engineer (e.g. in microwave 
cavities na ED, or is incompatible with the efficient 
operation of the metrological device (as in atom interfer¬ 
ometers HM3). 

In this paper, we present an alternative route to 
quantum-enhanced metrology based on purifications of 
mixed states. Physically, this involves entangling the 
probe with an auxiliary system before the probe is af¬ 
fected by </>, making measurements on both the probe and 
auxiliary system, and subsequently using correlations be¬ 
tween the two measurement outcomes in order to reduce 
the uncertainty in the estimated parameter (see Fig. [l]). 


This approach is advantageous in cases where it is easier 
to entangle the probe system with another system, rather 
than directly create highly entangled states of the probe 
system itself. An example of this is atom interferometry; 
although quantum squeezing can be produced in atomic 
systems via atomic interactions [T8H3T] , the technical re¬ 
quirements of high sensitivity, path separated atom inter¬ 
ferometers are better suited to enhancement via entan¬ 
glement with an optical system [321437] and information 
recycling 138H4T] . 

The structure of this paper is as follows. In Sec. [II 
we introduce in detail the central idea of this paper: that 
purifications of mixed states can possess a large quantum 
Fisher information (QFI), and therefore represent an ex¬ 
cellent resource for quantum metrology. In Sec. Ill we 


specialize to an TV-boson probe state and Mach-Zehnder 
(MZ) interferometer, and show how to engineer purifica¬ 
tions that yield sensitivities at and near the Heisenberg 
limit. Finally, in Sec. IV we present optimal measure¬ 
ment schemes that allow these quantum-enhanced sensi¬ 
tivities to be achieved in practice. 


II. QUANTUM FISHER INFORMATION FOR A 
PURIFICATION 

We can determine the best sensitivity possible for any 
given metrology scheme via the QFI, T 7 , which places 
an absolute lower bound on the sensitivity, A <f> >1 /, 
called the quantum Cramer-Rao bound (QCRB) [42H45] . 
This bound is independent of the choice of measurement 
and parameter estimation procedure, and depends only 
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FIG. 1. The unitary Uab — exp(— iHABt/h) entangles 
system A (probe) with system B (auxiliary) before system 
A passes through a measurement device described by TJ^ = 
exp(—z0GU). If measurements are restricted to system A, 
then the QFI for an estimate of 4> is Fa — F[Ga,Pa\, where 
pA = Tr^ {|Tas)(^as|}. If measurements on both sys¬ 
tems are permitted, then the QFI is Fab = F[Ga , |\1/ab)] = 
4Var {Ga)p a >F A - 


on the input state. Explicitly, if a state pA is input into 
a metrological device described by the unitary operator 
Ucf) = exp(— 1(\)Ga)i then the QFI is 


T a = T[Ga,Pa\ =2E ( a. (1 ) 

i,3 1 3 


where \ and |e^) are the eigenvalues and eigenvectors of 
Pa 5 respectively. If pA is pure, then Eq. 0 reduces to 
F a = 4Var(G A ). 

A naive consideration of the pure state QFI suggests 
that engineering input states with a large variance in 
Ga is an excellent strategy for achieving a high precision 
estimate of <fi. However, there are many operations on pA 
that increase Var(GU) at the expense of also decreasing 
the purity 7 = Tr {p^}. Since the QFI is convex in the 
state, any process that mixes the state typically decreases 
the QFI. Consequently, any improvement due to a larger 
Var(GU) is usually overwhelmed by reductions in the QFI 
due to mixing. 

In order to concretely demonstrate this point, we spe¬ 
cialize to an TV-boson state input into a MZ interferome¬ 
ter. As discussed in [46], this system is conveniently de¬ 
scribed by the SU( 2 ) Lie algebra [J*, Jj\ = ie^Jk , where 
eijk is the Levi-Civita symbol, for i = x,y, z. A MZ in¬ 
terferometer is characterized by Ga = J y , therefore, for 
pure states, a large QFI requires a large Var (J y ). 

Without loss of generality, we restrict ourselves to the 
class of input states 

Pa = dipV(A |a(§,<£>)}(a:(§, </>)|- (2) 

J 0 

Here |a(0, (p)) = exp(— i(pJ z ) exp (—iOJ y )\j,j) are spin co¬ 
herent states, where | j, m) are Dicke states with total 
angular momentum j = TV /2 and J z projection m. We 
focus on the following three states in class Q, which are 


Case (III) 



FIG. 2. (Color online) Husimi-Q function for Case (I) (a), 
Case (II) (b), and Case (III) (c). The projection in the J y 
basis, P(J y ) is shown for Case (I) (d), Case (II) (e), and 
Case (III) (f). N = 20 for all frames. 


in order of increasing Var (J y ): 
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Case (III): 

nv) = \ [<(*>-§)+<(*>+ 1 )]. 

(3c) 

These states can be conveniently visualized by plotting 
the Husimi-Q function [47l 138] 


Q(0, <p) = 2j 4 ^ 1 (a{0, y)\pA\a{9, p)) > ( 4 ) 


and the J y projection of the state, P(J y ) = (Jy\f>A\Jy), 
where J y \J y ) — J y \J y ) (see Fig. [ 2 ]) . 

None of these states yield sensitivities that surpass 
the standard quantum limit (SQL), A <j) = 1 /a/TV. In 
Case (I), pa is a pure spin coherent state, |o;( 7 r/ 2 ,0)), 
with Fa — 4Var (J y ) = TV. In Case (II), pA is an in¬ 
coherent mixture of Dicke states (i.e. it contains no off- 
diagonal terms in the | j,m) basis). Although 4Var(J y ) = 
TV (TV + l )/2 is much larger than for Case (I), the QFI 
is only Fa = TV/ 2 . Finally, Case (III) is an incoher¬ 
ent mixture of maximal and minimal J y eigenstates with 
4Var (J y ) = TV 2 , which is the maximum possible value 
in SU(2). However, since the state is mixed the QFI is 
significantly less than this, with Fa = TV/2. 

However, suppose the mixing in pA arises from en¬ 
tanglement with an auxiliary system B before system 
A passes through the metrological device (see Fig. [l]). 
Specifically, for an input pure state \^ab) of a compos¬ 
ite system A 0 T>, where pA = Tr# {\^ ab)^ ab\}, the 
QFI is 


Fab=P[G a ^ab)} 

= 4 {A!ab\G\\^ab) - {^abPa^ab) 2 ) 


= 4 ( Tta 
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Consequently, for a purification of pA the QFI only de¬ 
pends on the variance in Ga of pA HUES]. Our naive 
strategy of preparing a state with large Var(Gq) irrespec¬ 
tive of its purity is now an excellent approach. Indeed, 
in this situation the states in Cases (I)-(III) are now also 
arranged in order of increasing QFI, with Case (II) and 
Case (III) attaining a QFI of A(A + l)/2 and TV 2 , respec¬ 
tively. It is interesting to note that the QFI for Case (III) 
is the maximum allowable for A particles in SU(2) [50] . 
and is usually obtained via the difficult to generate spin- 
cat state, which is a macroscopic superposition, rather 
than a classical mixture, of spin coherent states. Note 
also that Tab is independent of any particular purifica¬ 
tion, and convexity implies that Tab > Aa. That is, in 
principle any purification of pA is capable of achieving 
sensitivities at least as good as, and usually much better 
than, pa itself. 

Quantum metrology with purifications is not simply 
a mathematical ‘trick’; physically, a purification corre¬ 
sponds to entangling the probe system A with some aux¬ 
iliary system B , and permitting measurements on both 
systems m- Therefore, the practical utility of our pro¬ 
posal depends crucially on the existence of an entangling 
Hamiltonian that can prepare pA in a state with large 
Var (GaW _ _ 

For the three cases described by Eq. ([2]) and Eqs. (|3]), 
a purification of pA can be written as 

3 

I’I'ab) = E c m\j, m) <g> \B m ), (6) 

m=—j 

with Case (I) corresponding to (B m \B n ) = 1, Case (II) 
corresponding to (B m \B n ) = £ n?m , and Case (III) corre¬ 
sponding to (B m \B n ) = 1(0) for \n — m\ even (odd). In 
the following section, we present a simple scheme that 
converts a shot-noise limited spin coherent state [such as 
Case (I)] to the enhanced QFI purifications of Cases (II) 
and (III). 


III. EXAMPLE ENTANGLING DYNAMICS 
LEADING TO INCREASED QFI 

Consider again the A-boson probe state (system A) 
input into a MZ interferometer (i.e. Ga — J y ). The QFI 


purification of pA can be written as 


T ab = 4 (<J y 2 

: ) - (Jy) 2 ) = To + A + 

(7) 

To = 

y(N + 2) — 2(jf), 

(8a) 

A = 

~ J-)) 2 , 

(8b) 

t 2 = 

-Al + J ->, 

(8c) 


0 ,1, and 2, in the J z basis; writing Tab in this form is 
very convenient for what follows. 

Before the interferometer, we assume the probe is cou¬ 
pled to some auxiliary system B via the Hamiltonian 

Bab = hgJzBs • (9) 

When system B is a photon field and Bb ^ proportional 
to the number of photons in the field, then Bab describes 
the weak probing of the population difference of an en¬ 
semble of two-level atoms with far-detuned light [371E1- 
ED, or dispersive coupling between a microwave cavity 
and a superconducting qubit [62li64] . We will explore 
this specific case shortly, however, for now we keep Bb 
completely general. If the initial system state is a prod¬ 
uct state |^ab(0 )) = |Tq) 0 after some evolution 

time the state of the system will be given by Eq. © with 
Cm = (m\^ A ) and \B m ) = exp(-imgt'H B )\'® b)- The 
reduced density operator of system A is then 

3 

Pa= X c n c* m C n - m \j,n)(j,m\, (10) 

n,m=—j 


where the coherence of system A is determined via 

Cn-m = (B m \B n ) = (g/sle-^-^WB |* B ). (n) 


When Cn-m = 1, the system remains separable and sys¬ 
tem A is a pure state, whereas if C n _ m = £ n?m then pA 
is an incoherent mixture of Dicke states. 

Using Eq. (10), To, T\, and T2 can be written as 


T 0 = f(N + 2)-2(J 2 ), 

(12a) 

T l =-(i(C l J + -ClJA)\ 

(12b) 

T 2 = -{C 2 (JI)+CUJ-)), 

(12c) 


where the above expectation values are calculated with 
respect to | Tq). The effect of the entanglement between 
systems A and B is entirely encoded in the coherences C\ 
and C2; coherences greater than 2nd order do not affect 
the QFI. 

Let us consider the effect on the QFI of each term in 
Eq. §. To is independent of the entanglement between 
systems A and B, and will be of order A 2 /2 if |Tq) has 
(A) rsj A (e.g. the spin coherent state |<a(7r/2, cp)) has 
(J 2 ) = A/4). This suggests that a sufficient condition 
for Heisenberg scaling is T\ ~ T2 ~ (J z ) ~ N. In fact, 
since T\ < 0, the maximum QFI state must necessarily 
have Ci = 0. In contrast, T 2 can be positive or negative, 
in which case a state with C 2 = 0 and another state with 
C 2 = 1 and T 2 ~ +A 2 /2 might both be capable of (near) 
Heisenberg-limited metrology. We consider examples of 
both states below. 


A. Case (II): Example dynamics yielding 

Tab ~ N 2 /2 


where J± = J x ± iJ y . Note that Ao, Jq, and T 2 depend To concretely illustrate the increased QFI a purifica- 
only on the matrix elements of pA with \n — m\ equal to tion of pA can provide, we assume system B is a single 
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bosonic mode, described by annihilation operator 6, and 
take Hb =b^b such that 

'Hab = hgJ z b^b . (13) 


If the initial state of system B is a Glauber coherent 
state \/3) [65], then the coherences described by Eq. (11) 
simplify to 


Cn—m — exp 


- \B\ 2 (l- e -i(n-m)gt^ 


(14) 


gt = 0 gt = 0.0251 gt = 0.0377 #£ = 0.15 






-50 0 50 -50 0 50 

Jy Jy 



-50 0 50 

Jy 


|C n _ m | 2 decays on a timescale gt ~ [(n — m) |/?| 2 ] -1 . Al¬ 
though the non-orthogonality of (/3\/3e' l °) ensures that 
C n _ m never actually reaches zero, it becomes very small 
for even modest values of \/3\ 2 . 

If the initial condition of system A is |4 /a) = | a(0, 0)), 
then Tab has the simple analytic form given by (see Ap¬ 
pendix [a| 

Tq — N (l+ ^ 2 ^ sin 2 o') ? (15a) 

T x = -N 2 sin 2 6 sin 2 (|/3| 2 sin (gt) + (f) e -^^^ 2 \ 

(15b) 

T 2 = N(1 2 N) sin 2 6 cos (|/3| 2 sin (2 gt) + 20) e ~ 2 ^ 2 sin ^ gt) . 

(15c) 


In contrast, calculating Ta via Eq. 0 requires the diag- 
onalization of pA , which must be performed numerically. 

We first demonstrate the effect of vanishing 1st and 
2nd order coherence on Tab by preparing system A in 
the maximal J x eigenstate, |a(7r/2, 0)), with N = 100, 
and a Glauber coherent state for system B with aver¬ 
age particle number \/3\ 2 = 500. The initial state for 
system A is precisely Case (I) [see Eq. (3a)], and has 
a QFI of N. As shown in Fig. [3j under the evolution 
of Eq. (Jl3|, pA tends towards an incoherent mixture of 
Dicke states [Case (II)], with the corresponding broaden¬ 
ing of the P(J y ) distribution. 


Figure [4][a) shows that both coherences C\ and C 2 
rapidly approach zero, which causes T\ and T 2 to van¬ 
ish [see Fig. [4](b)] . Consequently, Tab approaches To = 
N(N + l)/2, which allows a phase sensitivity of approxi¬ 
mately y/2 xHeisenberg limit [see Fig. [^c)]. In contrast, 
the effect of the mixing causes the QFI of pA itself to 
decrease from N to Ta = AT/2, with Ta < N for all t. 
This remains true even if the if Ga is rotated to lie in 
any arbitrary direction on the Bloch sphere. 

The oscillations in T\ and T 2 (and consequently Ta 
and Tab) before the plateau are due to the complex ro¬ 
tation of C\ and C 2 , which causes rotations of pA around 
the J z axis before being overwhelmed by the overall de¬ 
cay in magnitude. Furthermore, although the purity of 
the state also decays, it never vanishes, thereby illustrat¬ 
ing that it is not the entanglement per se that is causing 
the QFI enhancement for a purification of pA- 


FIG. 3. (Color online) Time snapshots of the Husimi-Q func¬ 
tion and J y projection illustrating the evolution of a maximal 
J x eigenstate under entangling Hamiltonian Eq. (13). The 
snapshots were chosen to correspond to times when the rota¬ 
tion around the J z axis is such that (J y ) = 0, which roughly 
corresponds to the local maxima of Tab in Fig. Be). (Pa- 
rameters: N = 100, \/3\ 2 — 500). 



gt 


FIG. 4. (Color online) Evolution of a maximal J x eigenstate 
under entangling Hamiltonian Eq. |l3] ). (a) Coherences |Ci | 2 
(blue solid line), IC 21 2 (red dashed line) and purity 7 (black 
dot-dashed line), (b): Three components of Tab [see Eq. ff\]: 
To (black dot-dashed line), T\ (blue solid line), and T 2 (red 
dashed line), (c) QFI for pa, Ta (blue dashed line), and a pu¬ 
rification of pa, Tab (red solid line). For reference, we have 
included N (black dotted line) and To = iV(iV + l )/2 ^ N 2 / 2 
(black dot-dashed line), which correspond to phase sensitivi¬ 
ties at the SQL and y/2x Heisenberg limit, respectively. (Pa¬ 
rameters: N = 100, \f3\ 2 — 500). 


B. Case (III): Example dynamics yielding 

Tab — N 2 


At gt = 7T, there is a revival in |C n | 2 for n even, but not 
for n odd. Figure [5] shows the Husimi-Q function under 
the evolution of Eq. for times close to gt = 7r, when 
the initial state of system A is the maximal J y eigenstate 
|a(7r/2,7r/2)). 

The QFI is initially zero, but the decay of C\ and C 2 
rapidly increases to Tab = To = N(N + l)/2 as in the 
previous example. As gt 7 r, the revival of IC 2 1 2 causes 
Tab to briefly increase to N 2 (see Fig. p|. This is the 
Heisenberg limit, which is the QFI of a (pure) spin-cat 
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gt = 0.977r gt = 0.9867T gt = 0.9927T gt = 7r 
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ing a Hamiltonian that does not conserve the J z projec¬ 
tion, such as when a spin flip in system A is correlated 
with the creation or annihilation of a quantum in system 
B. Here, we encounter scenarios where the interaction 
can either create or destroy coherences in the J z basis of 
system A, and although a significant QFI enhancement is 
still possible, it depends upon the initial state of system 
B. 

As a concrete illustration, consider the particle- 
exchange Hamiltonian 


FIG. 5. (Color online) Husimi-Q function and J y projection 
for an initial state |Tab( 0)) = |a(7r/2, 7 t/2 )) (g) \/3) under the 
evolution of Eq. ( fl3| for different values of gt. The Q function 
is symmetric about reflection of the J y axis, resulting in part 
of the function being hidden from view on the reverse side of 
the sphere. (Parameters: A — 100, \/3\ 2 = 500). 
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FIG. 6. (Color online) Evolution of a maximal J v eigenstate 
near gt — tt under entangling Hamiltonian Eq. (13). (a) \Ci\ 2 
(blue solid line), |C 21 2 (red dashed line), and Tfblack dot- 
dashed line), (b): Bo (black dot-dashed line), T\ (blue solid 
line), and T 2 (red dashed line), (c) Ta (blue dashed line) 
and Tab (red solid line). For comparison, we have included 
To — A(A +1)/2 N 2 /2 and the Heisenberg limit N 2 (black 
dot-dashed lines); Ta < A - for all t. (Parameters: A = 100, 
\/3\ 2 = 500). 


state and the maximum QFI for SU(2) [50]. At gt = 
7 r, pA is identical to a classical mixture of |a(7r/2, 7 t/ 2)) 
and |a(7r/2,—7 t/ 2)), however, its Q-function is similar to 
that of a spin-cat state, and purifications of it behave 
as a spin-cat state for metrological purposes. For these 
reasons, we call this state a pseudo-spin-cat state. 


C. Example dynamics for particle-exchange 
Hamiltonian 

In the previous two examples the J z projection was 
a conserved quantity, so any entanglement between sys¬ 
tems A and B can only degrade the coherence in the J z 
basis of system A (ultimately resulting in an enhanced 
QFI). The situation is more complicated when consider- 


n± = hg (j±6 f + j t 6) , 


(16) 


and assume that system A is initially prepared in the 
maximal J z eigenstate |4 /a ) = |a(0,0)) = | j,j) (N.B. this 
has Var(J y ) = A/4, and therefore a QFI of A). Then 
H- and 77+ physically correspond to Raman superradi¬ 
ance [551 EH] and quantum state transfer [34H371 M M 
processes, respectively. After some period of evolution, 
the combined state of systems A 0 B takes the form of 
Eq. (|6|. 

First, consider the case when the initial state for sys¬ 
tem A is a large amplitude coherent state (i.e. \^b) = 

| /?)). Here the addition/removal of a quantum to/from 
system B has a minimal effect on the state and the system 
remains approximately separable, since |(T> n |T> m )| 2 « 1. 
It is therefore reasonable to make the undepleted pump 
approximation 5—^/3, such that 77± —» hgf3J x (assuming 
/3 is real). Hence, the effect of the interaction is simply a 
rotation around the J x axis, which can create coherence 
in the J z basis, and so Ta = Fab < A for all time. 

In the opposite limit where the initial state of system 
B is a Fock state with A# particles, \^b) = |A#), then 


I B m ) = | A b ± (m - j)), 


(17) 


and (B m \B n ) = £ n?m . This ensures that the first and 
second order coherences vanish, and T\ — T 2 = 0 for all 
time. That is, as illustrated in Fig. [7| the state moves 
towards the equator and ultimately evolves to an inco¬ 
herent Dicke mixture [i.e. Case (II)]. As described in 
Sec. |III A[ and shown in Fig. [8j the QFI increases to 
a maximum of approximately Tab ~ A 2 /2. Although 
setting Nb = 0 (i.e. a vacuum state) leads to a larger 
variance in J 2 , Tab still reaches approximately 70% of 
A 2 /2. 


We therefore see that for the Hamiltonian (16), a large 


QFI enhancement is achieved provided the initial state 
| ^b) has small number fluctuations. Compare this to the 
Hamiltonian (13), where the choice \^b) — |A#) leads to 
no entanglement between systems A and T>, while in con¬ 
trast an initial state with small phase fluctuations (and 
therefore large number fluctuations), such as a coherent 
state, causes rapid decoherence in pA • 
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gt = 0 gt = 0.04 gt = 0.08 gt = 0.133 



FIG. 7. (Color online) Husimi-Q function and projection 
for an initial state |^as( 0)) = |a(0, 0)) 0 | Nb) under the 
evolution of 77- for different values of gt. (Parameters: N = 

100, N b = 20). 



FIG. 8. (Color online) Fa (blue dashed line) and Fab (red 
solid line) for an initial state |\I/ab( 0)) = |a(0, 0)} 0 | Nb) 
under the evolution of 77-. We have indicated N and TV 2 /2 
with black dotted lines for comparison. (Parameters: TV = 

100, N b = 20). 


IV. OPTIMAL MEASUREMENT SCHEMES 

Although the QFI determines the optimum sensitivity 
for a given initial state, it is silent on the question of 
how to achieve this optimum. It is therefore important 
to identify a) which measurements to make on each sys¬ 
tem and b) a method of combining the outcomes of these 
measurements - which we refer to as a measurement sig¬ 
nal (S) - that saturates the QCRB. We do this below for 
purifications of the incoherent Dicke mixture [Case (II)] 
and the pseudo-spin-cat state [Case (III)]. 


A. Optimal measurements for incoherent Dicke 
mixture [Case (II)]: 

It is worthwhile briefly recounting the optimal estima¬ 
tion procedure for a symmetric Dicke state | j, 0) input 
into a MZ interferometer. A MZ interferometer rotates 
J z according to J z {(j)) = U^J z U(j) = cos (j)J z — suuj)j x . 

Since symmetric Dicke states satisfy (J x ) = (J z ) = 
(J 2 ) = 0 and (J 2 ) = TV (TV + 2)/8, it is clear that the 
fluctuations in J z {(j)) contain the phase information, and 


therefore the quantity S = [ J z ( 0)] 2 oscillates between 
0 and TV(TV + 2)/8. It can be shown that at the op¬ 
erating point 0 —)> 0, Var(<S) 0, and the quantity 
(A (fl) 2 Var(<S)/(c^(<S)) 2 = 1 /.Fa, and therefore the 

signal S saturates the QCRB [H [67] [68:. 

For an incoherent Dicke mixture, we have (J x ) = 
(J y ) = (J z ) = 0, and (J 2 ) = TV(TV + l)/8. Unfortunately, 
the non-zero variance in J z (i.e. (J 2 ) = TV/4) implies 
that Var(<S) 0 for all 0, and the signal no longer sat¬ 
urates the QCRB. However, since the states | B m ) in the 
purification Eq. © are orthonormal, a projective mea¬ 
surement of some system B operator diagonal in the | B m ) 
basis projects system A into a J z eigenstate (i.e. a Dicke 
state). That is, these measurement outcomes on system 
B are correlated with J z measurement outcomes on sys¬ 
tem A. Therefore, subtracting both measurements yields 
a quantity with very little quantum noise. 

More precisely, if we can construct an operator Sb on 
system B that is correlated with J z measurements on 
system A (i.e. = J z \^ab}), then we can con¬ 

struct the quantity So = J z — Sb which has the property 
(<So) = (Sq) = 0. This motivates the signal choice 

S = (u^SoU^j = (cos (f)J z - Sin <j>j x - S B ) 2 - (18) 

Using SbI'&ab) = J z \^ab) and the fact that non-J^ con¬ 
serving terms vanish due to the absence of off-diagonal 
terms in the J z representation of (e.g. expectation 
values with an odd power of J x vanish), we can show that 


(S) = (J 2 ) (cos 0 — l) 2 + (J 2 ) sin 2 0, (19a) 

(^ 2 ) = (Jz )( cos 0 “ !) 4 + (Jx) sin4 ^ 

+ ( J Z J X + J X J Z + 4 J Z J X J X J Z ) sin 2 0(cos (j) — l) 2 
+ 2 i((J z J X Jy — JyJxJz )) sin 2 (p COS (p (cos — 1) 

+ ( Jy) cos 2 0sin 2 (p. (19b) 

Note that the above expectation values can be taken with 
respect to | ^ab) o r pA- The best sensitivity occurs at 
small displacements around 0 = 0. Taking the limit as 
(f) — y 0 and noting that (J 2 ) = (J 2 ) gives 


(A0) 2 = 


Var(<$) 

(^(<S » 2 


( f >=0 




1 

Fab 


( 20 ) 

This demonstrates that the signal Eq. (18) is optimal 
since it saturates the QCRB. 


B. Optimal measurements for pseudo-spin-cat 
state [Case (III)]: 

Pure spin-cat states have the maximum QFI possible 
for TV particles in SU(2), are eigenstates of the parity 
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operator, and indeed parity measurements saturate the 
QCRB |69|. Pseudo-spin-cat states (case (III)) also have 
maximal QFI, and since (B n \B m ) = 1(0) for \n — m\ 
even(odd), a projective measurement of system B yields 
no information other than the parity of the J z projec¬ 
tion. This suggests that a measurement of parity could 
be optimal. 

In analogy with Case (II), our aim is to construct an 
operator So where the correlations between systems A 
and B lead to a reduction in Var(<S 0 ) and the system 
mimics a pure spin-cat state. Introducing the quantity 

<5o = UaSb IIaIIb, (21) 


where ILa(b) is the parity operator for system A(B ), 
defined by U A \j,m) = (-l) m | j,m) and IL B \B m ) = 
(—l) m |H m ), we see that pseudo-spin-cat states satisfy 
= |^ab), and therefore Var(<So) = 0. This mo¬ 
tivates the signal choice S = U^SoU#- 

To calculate the sensitivity, we need to compute (S) 
and (<S 2 ). Trivially, (<S 2 ) = 1 for all states. For (/> 1, 

expanding to second order in (j> gives 


(5) * ((1 + i^Jy ~ WJ l y )So(\ - i^Jy ~ \VJ 2 )) 
= 1 icj) (^(JySo) — (SoJy)^ 

(JySoJy) - \ (< J 2 y S 0 ) + (S 0 J 2 )) 


+ <f 
+ o(A), 


( 22 ) 


The relation (B n \B n ±i) = 0 ensures that terms linear in 
J y go to zero: 


(J+) ^ ^ c m c n \ji m) {Bn \Bm) 

m, n 

OC E C n ^n,m+1 (B n \B 

m) 

n, m 

= c m c m+ i {B m + 11 B m ) = 0. (23) 

rri 


However, unlike Case (II), the condition ( B n \B n ± 2 ) = 1 
preserves terms such as {J+). Noting that J y flips the 
parity of any state in subsystem A but not subsystem B: 


^aJ V \^Ab) = —Jy^A^AB), (24a) 

A B Jy\^ Ab) = Jyt^B^AB), (24b) 

and using 5 0 |^ab) = gives 

(JySoJy) = ~(SoJ 2) = -(JySo) = ~(J 2 y). (25) 


Since Sq = 1 implies that S 2 = 1, we obtain 

Var(<S)=40 2 (J 2 ) + (9(</> 4 ), (27) 

and consequently 

(a , )2 = Var (5) = 4^> 2 (J 2 ) = _J_ = J_ 

mS)) 2 16 </> 2 (J 2 ) 2 4(J2) Tab' 

(28) 

This demonstrates that the signal saturates the QCRB 
and is therefore optimal. 

The optimal estimation schemes presented in 
Secs. uva] and hvh illustrate a somewhat counter¬ 
intuitive fact: although the optimal measurement of 
system B for a pseudo-spin-cat state provides less 
information about system A than for an incoherent 
Dicke mixture, the pseudo-spin-cat state yields the 
better (in fact best ) sensitivity. 


C. System B observables that approximate optimal 
measurements 

We now turn to the explicit construction of physical 
observables that approximate In general, the choice 
of Sb depends upon the specific purification of pA • Phys¬ 
ically, the initial state of system B and the entangling 
Hamiltonian matter. However, there is no guarantee that 
Sb exists, and if it does there is no guarantee that a 
measurement of this observable can be made in practice. 
Nevertheless, as we show below, it may be possible to 
make a measurement of an observable that approximates 
Sb , and can therefore give near-optimal sensitivities. 


1. Case (II) 


To begin, consider the situation in Sec. Ill A the evo¬ 
lution of the state |a(7r/2, 0)) ® \/3) under the Hamilto¬ 
nian (13). We require Sb\^ab) = JzI'&ab)- After some 
evolution time t : 


\9ab)~ E c m \j,m)\Pe~ im9t ). (29) 

m=—j 

Clearly, the phase of the coherent state is correlated with 
the J z projection of system B. This can be extracted 
via a homodyne measurement of the phase quadrature 
Yb = i(b — id) [70]. In fact, provided mgt <C 1, phase 
quadrature measurements of \/3 exp(— imgt)) are linearly 
proportional to the J z projection: 

((3e- irngt \Y B \(3e- irngt ) = 2/3sin (mgt) « 2 (3mgt, (30) 

where without loss of generality we have taken /3 to be 
real and positive. Consequently, the scaled phase quadra¬ 
ture 


Sb 


Y b 

2f3gt 


Therefore 


{S) = l-2<p 2 {j 2 y )+0{4?). 


(26) 


(31) 















satisfies 


(^ab | (j z ~ Sb ^ | ^ab) ~ 0, (32a) 

(^Ab| (tz - 5 b) l^yis) ~ / 2j gfft)2 ’ ( 32b ) 


and so the fluctuations in (J z — Sb) become arbitrarily 
small (and <S# becomes perfectly correlated with J z ) as 
((3gt) 2 becomes large. This suggests that Eqs. (18) and 
( |3T| ) should be a good approximation to an optimal mea¬ 
surement signal. 

More precisely, assume that 



< f> (rads) 


1 Ngt 

<<l 


(33) 


The first inequality ensures that (/3e~ in9t \l3e~ irn9t ) « 
d n?m and so is approximately an incoherent Dicke mix¬ 
ture, while the second inequality implies that we are in 
the linearized regime where Eq. (30) and Eqs. (32) hold. 
Then the signal 


^approx — {jz COS 0 Jx sin 0 2,gt 


(34) 


yields the sensitivity 

1 


(A <Pf 


+ 


(W» 

4 

(2 Pgt) 2 


|Var(<S) 


sin 2 (0/2) 


At) 


P 2 x z ' (2 pgty 

(cos<(>-1) 2 (J 2 )+sin 2 0(j|) ), 

(35) 


where Var(<S) and 8^(3) are given by the expecta¬ 
tions p9| of the optimal signal Eq. (18). Figure ([9| shows 
Eq. (|35|) compared to an exact numeric calculation. 

Condition ( [33] ) typically ensures that the term propor¬ 
tional to (J 4 ) is small in comparison to the term propor¬ 
tional to 1/(2 ftgt) 2 . We therefore see that our approx¬ 
imate signal <S ap prox gives a sensitivity worse than the 
QCRB, and furthermore at an operating point 0^0. 
Nevertheless, A0 approaches the QCRB at 0 = 0 as 
/ 3 2 and (2 ftgt) 2 approach infinity. Therefore, for a suf¬ 
ficiently large /3gt, we can achieve near-optimal sensi¬ 
tivities close to (j) = 0. This is illustrated in Fig. [9] 
When /3gt = 10, we find that A <fi is very close to the 
QCRB. In contrast, for /3gt = 1, the imperfect corre¬ 
lations between J z and Sb prevent the sensitivity from 
reaching the QCRB; nevertheless, the sensitivity is still 
below the SQL. Note that there is a slight deviation be¬ 
tween Eq. (35) and the numerical calculation of the sensi¬ 


tivity using the state (29). This is due to terms neglected 


by our approximations; in particular, the nonlinear terms 
ignored by linearizations such as Eq. ( |3()| ) and those ne¬ 
glected terms that arise due to the small (but strictly 
non-zero) off-diagonal elements of pA • 


FIG. 9. (Color online) A0 versus 4> using Eq. ( [34] ) for a state 
of the form Eq. ( [ 29 ] ) , with N — 100, and gt = 10 -2 . The 
blue dot-dashed line is with \/3\ 2 — 10 6 (/ 3gt = 10), and 
the red solid line is for \/3 \ 2 = 10 4 (/ 3gt = 1). The red 
dashed line shows the approximate expression for the sen¬ 
sitivity (Eq. (35)) for \/3 \ 2 = 10 4 . For \/3 \ 2 — 10 6 , the nu¬ 
merical calculation and Eq. (35) are identical. The upper and 
lower black dotted lines represent the standard quantum limit 
(1/y/N), and y/2/N respectively. The divergence in A <f> close 
to 0 = 0 in both cases is due to the imperfect correlations 
between Sb and J z leading to non-zero variance in S. If the 
correlations were perfect and Var(<S)|$—o = 0, A (j) would reach 
exactly 1 j \] Tab at <f> = 0. 


2. Case (III) 


Now, consider the situation in Sec. mm the evolution 
of the state |a(0, 0)) ( 8 ) \/3) under the Hamiltonian (13) 
that at gt = 7r approximately results in a pseudo-spin- 
cat state. 

In order to find an operator that approximates Sb = 
II£, we introduce the amplitude quadrature operator 
Xb = (5 + m), and notice that 


(Pe- im *\X B \Pe 


r ) = 2/3(—l) m 

= 2P(j,m\il A \j,m). 


(36) 


That is, amplitude quadrature measurements of \/3e~ irn7T ) 
are proportional to parity measurements on system E>, 
which are directly correlated with parity measurements 
on J z eigenstates. Indeed, the quantity 


S 0 = il A ^, (37) 

has a variance Var(<So) = 1/(2 /3) 2 that becomes vanish¬ 
ingly small as the amplitude of the coherent state is in¬ 
creased. We therefore expect the signal 

Approx = U\ y (38) 

will be a good approximation to the optimal measure¬ 
ment S. 
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FIG. 10. (Color online) Phase sensitivity of the approximate 
signal Eq. ( [38] ) for a state of the form Eq. ( [29] ) at gt — tt (i.e. 
a pseudo-spin-cat state) with N = 20. The blue solid line and 
red dashed line are for \/3\ 2 = 30 and \/3\ 2 = 5, respectively. 
The black dotted line indicates the Heisenberg limit Acj) — 
1/N (which is the QCRB). Note that the vertical axis is a 
linear scale. 


Fi gure 10 shows the sensitivity for a state of the form 
Eq. (29) at gt = tt with N = 20. When \/3\ 2 = 30, the 
sensitivity is very close to the Heisenberg limit, while 
for |/31 2 = 5 there is a slight degradation in the sensi¬ 
tivity due to imperfect correlations. In contrast to the 
approximate optimal measurement scheme for Case (II), 
which requires a large amplitude coherent state, here the 
signal (38) is almost optimal even for small amplitude co¬ 
herent states. This is because (/3e~ l7r \/3) = exp(—21/51 2 ) 
is approximately zero even for modest values of /3. 

In situations where system A is an ensemble of atoms, 
and system B is an optical mode, it would be challenging 
to achieve the strong atom-light coupling regime required 
for gt = it. On the other hand, the choice of an initial 
coherent state for system B ensures that the sensitivity 
is reasonably insensitive to losses in system B. In partic¬ 
ular, since particle loss from a coherent state acts only to 
reduce the state’s amplitude, provided the coherent state 
remains sufficiently large after losses in order to satisfy 
the requirements for near-optimal measurements, near- 
Heisenberg-limited sensitivities should be obtainable. 


3. Particle-exchange Hamiltonian 


Finally, for completeness we include the optimal mea¬ 
surement scheme for the state attained after evolving the 
product state | q(Q, 0))<S>\Nb) under the Hamiltonian (16) 
(see Sec. Ill C|) . The optimal measurement signal is sim¬ 
ply Eq. (18) with 


Sb = y ± (iV b - &*&) . (39) 

This choice of Sb can be constructed by counting 
the number of particles in system 13, and it satisfies 
Sb\^ab) = JzY&ab) as required. As any entangling 


Hamiltonian of the form 

n± = (j ± St + j T bf + Y / B j ,kJz (^y (40) 

k j,k 

will lead to a state of the form Eq. with | B m ) given 
by Eq. ( [T7| ) (assuming an initial state \S&a) — |ck(0,0)), 

b) = pYe)), Eq. ( [39] ) also transfers to these systems. 

V. DISCUSSION AND CONCLUSIONS 

We have shown that purifications of mixed states rep¬ 
resent an excellent resource for quantum metrology. In 
particular, we showed that if probe system A and auxil¬ 
iary system B are entangled such that the 1st and 2nd or¬ 
der coherences of system A vanish, then near-Heisenberg- 
limited sensitivities can be achieved provided measure¬ 
ments on both systems A and B are allowed. Although 
we focused on the situation where this entanglement is 
generated via a few specific Hamiltonians, our conclu¬ 
sions hold irrespective of the specific entanglement gen¬ 
eration scheme. 

While preparing this paper, we also numerically exam¬ 
ined the effect of decoherence on the sensitivity of our 
metrological schemes. In particular, we found that the 
effect of particle loss, spin flips, and phase diffusion on 
purifications of the pseudo-spin-cat state from Fig. [2] was 
identical to that of a pure spin-cat state. 

Although these purified states are no more or less ro¬ 
bust to decoherence than other nonclassical pure states, 
there are situations where they are easier to generate. 
The example we are most familiar with is atom interfer¬ 
ometry, where atom-light entanglement and information 
recycling is more compatible with the requirements of 
high precision atom interferometry than the preparation 
of nonclassical atomic states via interatomic interactions 
[39] . However, controlled interactions are routinely en¬ 
gineered between atoms and light [37 JTTH73] . supercon¬ 
ducting circuits and microwaves [71175], light and me¬ 
chanical systems G5I, and ions and light IZZHZ3I. Given 
that high efficiency detection is available in all these sys¬ 
tems [801485] . the application of our proposal to a range 
of metrological platforms is plausible in the near term. 

It is important to note that although the QFI ap¬ 
proaches the Heisenberg limit (Fab = 3V 2 , in Case (III), 
for example), this is not the true Heisenberg limit, as 
N refers only to the number of particles in system A 
(which pass through the interferometer), rather than the 
total number of particles N t in system A and system B. 
However, there are some situations where the number of 
particles in system A is by far the more valuable resource, 
which is why it makes sense to report the QFI in terms of 
N rather than N t . For example, consider the case of iner¬ 
tial sensing with atom interferometry, where system A is 
atoms, and system B is photons. The atoms are sensitive 
to inertial phase shift, but it is difficult to arbitrarily in¬ 
crease the atomic flux. However, a gain can be achieved 
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by adding some number of photons to the system, which 
are comparatively ‘cheep’ compared to atoms. 

Finally, we note that not all quantum systems are cre¬ 
ated equal; certain quantum information protocols, such 
as quantum error correction [86j and no-knowledge feed¬ 
back m, are better suited to some platforms than oth¬ 
ers. Our proposal allows an experimenter to both per¬ 
form quantum-enhanced metrology and take advantage 
of any additional benefits a hybrid quantum system pro¬ 
vides. 
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Appendix A: Derivation of Eqs. © 

Here we derive the QFI Tab for a MZ interferometer 
with the following entangled input: 

\*AB(t))=e- i ^»\d,A®\P), (Al) 

where Ny = tfb, system B is initially in a coherent state 
| /3), and system A is initially in a spin coherent state 
1 0,p). Any spin coherent state can be defined by rotating 
the maximal Dicke state on the top pole of the Bloch 
sphere an angle 0 about the J y axis and an angle p about 
the J z axis: 

\0,p) = R(0,A\j,j) = j,j). (A2) 

Recall that j = N/ 2, where N is the total number of 
system A particles. 

The QFI is 

Tab = 4Var A'j , (A3) 

where the expectations in the variance are taken with 
respect to the initial separable state, 1 0,p) ( 8 ) \/3). 

By application of the Baker-Campbell-Hausdorff for¬ 
mula 


1 r 


e A Be- A = B+[A,B] + -\A, [.A,B} 




A [A,B] 


+ 


(A4) 


Therefore, since the initial state is separable, we obtain 
Tab = (sin 2 (gtN b )) (J 2 X ) + (cos 2 {gtN b )) (J 2 V ) 

+ (sin (gtNfr) COs(gtNy)) (J x Jy J - JyJx) 

- ({sm(gtN b ))(J x ) + (cos(gtN b )){Jy)^j . (A6) 

The system A expectations are more easily computed 
by rotating the operators by R(6 , p) and then taking ex¬ 
pectations with respect to the Dicke state | j,j). Specifi¬ 
cally, by virtue of 

W{e,<p)j x R[e,ip) = cos 0 cos (pJ x + cos 0 sin pj y 



— sm6J z , 

(A7a) 

w(0, p)J y R(6 , (p) = — sin^J^ + cos pJ y , 

(A7b) 

and the application of J± = J x db iJ y with 

J±\j,m) = y/j{j + 1 ) - (m±l)\j,m±l), 

(AS) 

we obtain 

(j x ) = j sin 0 cos p, 

(A9a) 


(J y ) = j sinOsinp , 

(A9b) 


(Jx) = ^ (! + ( 2 i - !) sin 2 9 cos 2 p) , 

(A9c) 


(Jy) ~ ^ (l + (2j 1) sin 2 9 sin 2 p) , 

(A9d) 

( J x J y + J y J x ) = j (2 j — 1 ) sin 2 0 sin p cos p. 

(A9e) 

With some 

simplification this gives 



Fab = 2j(l + sin 2 6[(2j - l)(sin 2 (gtN b + p)) 

~ 2 j (sin (gtN b + p)) 2 ] j. (A10) 

Incidentally, by setting t = 0 we can see that the QFI for 
a spin coherent state input never exceeds the standard 
quantum limit: 

T[J y , |0, p)] = 2 j (l — sin 2 0 sin 2 p) < N. (All) 

In order to compute the system B expectations, note 
that 

(sin (gtN b + p)) = -~ ((e^Ro+v)) _ (c^W+f))) 

(Al2a) 

2 2 — _ ^g-2 i(gtN b +<p)^ 

4 (Al2b) 


Furthermore, for any m, 

^ e im(gtN b +(p) ) _ ^nrup 

= exp [-|/?| 2 (1 - e im9t ) + 2 mp\ . (A13) 

Substituting Eqs. ( A12| ) and Eq. ( |A13[ ) into Eq. (A10) 
gives 


it can be shown that 

e igtJ z N b j ye -igtJ z N b = S m(gtN b )J x + cos{gtN b )J y . (A5) 


Fab = Fq + F-\ + F 2 , (A14) 


with the expressions for To, 
Eqs. ([15). 


T \, and T 2 listed in 
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